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ABSTRACT 

We examine large-scale fluctuations in the He II Lya forest transmission during and after He II 
rcionization. We use a simple Monte Carlo model to distribute quasars throughout a large volume 
and compute the resulting radiation field along one-dimensional skewers. In agreement with previous 
studies, we find that the rarity of these sources induces order unity fluctuations in the mean optical 
depth after reionization, even when averaged over large segments (~ 10-100 Mpc across). We compare 
our models to existing data along five He II Lya forest lines of sight spanning z ~ 2-3.2. The large 
cosmic variance contained in our model plausibly explains many of the observed fluctuations at z < 2.5. 
But our models cannot accommodate the large fluctuations toward high optical depths on ~ 30 Mpc 
scales observed at z ~ 2.7-2.9, and the measured optical depths (r c ff > 4) at z > 2.9 are difficult 
to explain with a smoothly-evolving mean radiation field. In order to better understand this data, 
we construct a toy model of He II reionization, in which we assume that regions with the smallest 
radiation fields in a post-reionization Universe (or farthest from strong ionizing sources) are completely 
dark during reionization. The observed fluctuations fit much more comfortably into this model, and 
we therefore argue that, according to present data, He II reionization does not complete until z < 2.9. 
Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 

Our premier tool for studying the intergalactic medium 
(IGM) is the Lya forest, from which we can measure the 
absorption properties of the low-density "cosmic web" 
lying between galaxies. Although observations of the 
H I forest have a long and rich forty-year history (Rauch 
1998), recently there has been great progress in mea- 
suring absorption from other ions, including both heavy 
elements and helium, the second most common element 
in the Universe. 

Unfortunately, observing the He II Lya forest is dif- 
ficult, because the transition's rest wavelength of 304 
A places much of the absorption in the far ultraviolet, 
even for material at z ~ 2-3. Not only do ultravi- 
olet telescopes present formidable technical challenges, 
but studying the forest also requires a luminous far-UV 
source to provide a background light. Such lines of sight 
are rare, because any that intersect strong H I absorbers 
do not retain their UV flux. To date, five quasars at 
z ~ 2-3.5 have satisfied these criteria and have high- 
resolution He II Lya forest measurements (see §5 below 
for references to the individual quasars), although very 
recently more than twenty other potential targets have 
been identified (Syphers et al. 2009a, b). 

Nevertheless, the known lines of sight are particu- 
larly interesting, because current data (from these mea- 
surements and other indirect techniques) suggest that 
quasars fully ionized He II by z ~ 3. Most obviously, 
there is a substantial decrease in the mean Lya forest 
absorption beyond z ~ 2.9 (Dixon & Furlanetto 2009), 
which is compatible with models where large swathes of 
He II are cleared out by quasar light at about that time 
(Sokasian et al. 2002; Furlanetto & Oh 2008b; McQuinn 
et al. 2009). Indirect evidence from the IGM temper- 
ature also suggests reionization at z ~ 3.2-3.4 (Schaye 
et al. 2000; Ricotti et al. 2000; Lidz et al. 2009), which is 
consistent with models of helium reionization ending at 



these slightly higher redshifts (Gleser et al. 2005; Furlan- 
etto k Oh 2008a; McQuinn et al. 2009). 

Here, we will study the variability in the He II Lya 
forest data and its consequences for reionization. Specif- 
ically, recent work has shown that the high-energy radia- 
tion background fluctuates strongly both during and af- 
ter helium reionization, thanks largely to the rarity of the 
bright quasars providing the ionizing photons (see Fardal 
et al. 1998; Maselli & Ferrara 2005; Bolton et al. 2006; 
Meiksin 2007; Tittley & Meiksin 2007; Furlanetto 2009a 
for theoretical work and Shull et al. 2004; Fechner et al. 
2006; Fechner & Reimers 2007a for observations). More- 
over, these fluctuations can span large physical scales 
(Furlanetto 2009b). We might therefore expect that cos- 
mic variance from the small number of well-studied lines 
of sight may limit any firm conclusion based on the op- 
tical depth. 

On the other hand, the existing observations already 
exhibit substantial fluctuations and may themselves sug- 
gest something useful about reionization. We naturally 
expect larger fluctuations during reionization, because 
the radiation field can be extremely large near quasars 
but almost zero in the He II regions far from these sources 
(Meiksin 2007; Furlanetto 2009a). Here, we will quan- 
tify how these fluctuations impact measurements of the 
mean Lya forest absorption and show that the observed 
variations are by themselves powerful probes of the radi- 
ation field. We will use a simple Monte Carlo model and 
ask whether the existing data are consistent with He II 
reionization ending at z ~ 3. Our approach is similar 
to Fardal et al. (1998), although we use updated input 
parameters, add a toy model for reionization itself, and 
have the advantage of comparing to significantly more 
data. 

This paper is organized as follows. In §2, we de- 
scribe the Monte Carlo model that we use for the post- 
reionization limit. We present a simple toy model of He II 
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rcionization in §3 and illustrate the basic characteristics 
of these models in §4. In §5, we compare our results to 
existing data. Finally, we conclude in §6. 

In our numerical calculations, we assume a cosmol- 
ogy with O ro = 0.26, Q A = 0.74, tt b = 0.044, H = 
100ft, km s- 1 Mpc" 1 (with h = 0.74), n = 0.95, and 
£78 = 0.8, consistent with the most recent measurements 
(Dunkley et al. 2009). Unless otherwise specified, we use 
comoving units for all distances. 

2. MONTE CARLO MODEL 

2.1. The Post-Reionization Limit 

We have previously examined the distribution of ioniz- 
ing background intensities in the post-reionization limit 
with an analytic model (Furlanetto 2009a; see also Zuo 
1992a; Meiksin & White 2003; Meiksin 2007). This 
model predicts substantial line-of-sight variations in the 
specific intensity J, even on large scales, which can be 
quantified by, for example, the power spectrum of the 
ionizing intensity (Zuo 1992b; Furlanetto 2009b). How- 
ever, the power spectrum does not perfectly match the 
method used to analyze the observations, where optical 
depth measurements are binned over large scales in or- 
der to reduce the (substantial) errors in measurements 
of individual features. These differences include: (1) 
real quasar absorption spectra allow only a finite dy- 
namic range to measure transmission, (2) observations 
are typically binned with sharp filters in real space, (3) 
the (measurable) ionization rate requires an integral over 
the entire quasar spectrum, rather than evaluation at a 
single frequency, and finally (4) observations are sensitive 
to the density field and temperature of the IGM in addi- 
tion to the radiation background studied by Furlanetto 
(2009b). 

To address (most of) these concerns, we instead use a 
simple Monte Carlo model to construct the ionizing field 
along example lines of sight. Our model is very similar to 
Fardal et al. (1998), except that we use more recent data 
as inputs to the calculation. We begin by populating a 
box of size ~ 1 Gpc 3 with quasars by randomly sampling 
the Hopkins et al. (2007) luminosity function at a given 
rcdshift; the precise size of the box varies with the appli- 
cation but is always at least this large. 1 We also draw the 
far-UV spectral index of each quasar (assumed to be a 
power law at A < 1050 A) from a Gaussian distribution 
with (a) = 1.5 and variance a 2 — 1, except truncated 
between a = (0.5, 3.5). This is roughly consistent with 
the observed distribution of Telfer et al. (2002), and the 
central value is also consistent with the recent direct mea- 
surements of Syphers et al. (2009b). We assign quasar 
positions randomly throughout the box. Note that this 
procedure does not account for quasar clustering, which 
is substantial (Shen et al. 2007) and will tend to amplify 
the fluctuations (see discussion in §6). 

Next, we choose a line of sight of a specified length. To 
do this, we randomly choose the segment's starting posi- 
tion within the box and its direction on the unit sphere, 
wrapping the line of sight around the box if necessary. 
We then compute the ionizing intensity in 1 Mpc steps 
along the line of sight. In each case, we sum over all of 

1 To generate all of our random numbers, we use the Mcrscnnc 
Twister algorithm, which has a period of at least 2 19937 — 1 (Mat- 
sumoto & Nishimura 1998). 



the quasars in the box, 2 so that the specific intensity at 
a given frequency (in units of erg cm~ 2 s _1 Hz -1 ) is 
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where J* = L^/ (4ttRq), Li is the quasar's luminosity, 
L^ is the mean luminosity, is the distance between 
the i th quasar and the point of interest, and Rq is the 
attenuation length of the ionizing photons (see below). 

Finally, we are generally interested in the ionization 
rate T rather than the specific intensity at a single fre- 
quency, 



r = / dv J v an e n{y), 
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where (THoii is the photoionization cross-section. If we 
approximate <7H e ii oc v~ z relatively close to the ioniza- 
tion threshold (Verner et al. 1996) and consider a quasar 
with spectral index at, we have 
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where r, is the ionization rate from the ith quasar, the 
last factor accounts for the spectral decline between 1050 
A and the He II ionization edge, and T denotes an aver- 
age over the quasar distribution. Note that here we have 
assumed that Rq is independent of frequency. In reality, 
the attenuation length can depend strongly on frequency: 
for example, Rq cx l/a cx v z in a uniform medium. The 
dumpiness of the IGM moderates this somewhat; for ex- 
ample, a column density distribution cx implies 
Rq cx v 15 (Paresce et al. 1980). Our simple treatment 
is only valid in the limit that all the absorption comes 
from extremely high-column density systems. However, 
a more realistic approach will only decrease the fluctu- 
ations and hence strengthen our conclusions, because it 
would increase the attenuation length at large frequencies 
and make the observed fluctuations even more difficult 
to explain. 

For our calculations here, we generate three indepen- 
dent boxes at each relevant redshift. Within each box, 
we generate 100 lines of sight, each ~ 1 Gpc in length. 
This procedure provides 3,000 samples of ~ 100 Mpc 
segments (or many more for comparisons to shorter seg- 
ments). We choose the precise dimensions to match the 
data points in the observed sample. In practice, there is 
relatively little difference in the results between the three 
independent boxes. 

The bottom panel of Figure 1 shows three example 
lines of sight generated with this method; here, we have 
assumed z — 3 and R = 30 Mpc (see below) . We show 
the ionization rate scaled to its mean value, (r) = 3-/VrjI\, 
where No is the mean number of quasars within one at- 
tenuation volume and is the ionization rate for a point 
with specific intensity J* (Meiksin & White 2003). These 
examples illustrate several important aspects of the sig- 
nal (see also Furlanetto 2009b): (1) the correlations can 
extend over extremely large scales, as in the long trough 

2 In detail, we use the closest position of each quasar in a periodic 
tiling of the box, so that we effectively include all sources within 
1 Gpc of every point. This scale is much larger than the assumed 
attenuation length, so more distant sources can safely be ignored. 
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Fig. 1. — Ionization rate (scaled to its mean) along three example 
lines of sight from our Monte Carlo simulations. The bottom panel 
shows the post reionization limit; the top panel shows the same 
lines of sight processed with our toy model of He II reionization, 
with 20% of the pixels assumed to be outside of illuminated He III 
regions. All take Rq = 30 Mpc. 



at ~ 200-500 Mpc in the dashed curve; (2) rare small- 
scale "fluctuations" are present, because of the 1/r 2 pro- 
files around individual sources (these dominate the power 
spectrum); (3) the median of J is a few times smaller 
than the mean; and (4) fluctuations far above the me- 
dian level are relatively common, but the ionization rate 
never falls significantly below that level. This last point 
is because the accumulated contributions of the many 
quasars at distances larger than R Q provides an effective 
floor to r, and it is crucial for the following. 

2.2. From the Ionizing Background to the Effective 
Optical Depth 

We now wish to transform our results to the observable 
He II Lya optical depth. Because we are interested in 
spatial fluctuations over relatively large scales, we will 
average the observed transmission T over large segments, 
parameterizing it in terms of the effective optical depth 

7"off , 

T cff = -lnX. (4) 

Here T denotes the average transmission over a single 
contiguous segment. We will use to (t c s) to denote the 
effective optical depth when averaged over many inde- 
pendent segments. 

We therefore need a prescription relating T to the local 
forest transmission. If the IGM were a uniform medium, 
the relation would be trivial: in ionization equilibrium, 
the He II fraction obeys XHcii oc r _1 , so t a T" 1 as 
well. However, in the real IGM a forest of discrete clouds 
provides most of the absorption. In this case, 
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where d 2 N / dN^ciidz is the number of absorbers per 
unit redshift and per unit column density and W\ = 
f(l — e~ T )d\ is the equivalent width of each line. Be- 
cause the abundances of lines in the He II forest are 



poorly constrained at present, we write iVHeii = V^m 
and use the observed column density distribution of the 
H I forest to compute the integral; in the optically thin 
regime (relevant for the vast majority of low-column 
density systems providing Lya absorption), rj is inde- 
pendent of -/VhcII (Haardt & Madau 1996; Fardal et al. 
1998; Faucher-Giguere et al. 2009). For such systems, a 
good approximation to the column density distribution 
is oc N^f , where f3 ~ 1.5 (Fardal et al. 1998, and refer- 
ences therein). In that case, the integral in equation (5) 
can be rewritten as (Fardal & Shull 1993), 
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where (r) is the mean ionization rate. Note that j3 = 
can mimic a uniform medium. In fact, existing obser- 
vations of the He II forest are not of sufficient signal-to- 
noise to rule out the presence of a uniform component in 
the overall absorption, although most of the absorption 
does clearly come from discrete systems that can be as- 
sociated with corresponding H I absorbers (e.g., Fechner 
et al. 2006). Fardal et al. (1998) showed that a com- 
pletely uniform IGM roughly doubles fluctuations in the 
optical depth. 

Of course, T does vary spatially throughout the simu- 
lation. We, therefore, compute the transmission (relative 
to the mean) separately in each 1 Mpc pixel of the line 
of sight. We then set the proportionality constant in 
equation (6) so that the mean transmission (r e ff) has a 
specified value when averaged over our entire sample of 
lines of sight (usually chosen to match observations or 
an input model). Finally, we compute t c s averaged over 
larger continuous distances (typically ~ 30-100 Mpc) by 
averaging over all the pixels within each such region. 

2.3. Other Sources of Fluctuations 

In the previous subsection, we showed how T induces 
fluctuations in the transmission. However, the He II Lya 
forest is sensitive not only to the local radiation field but 
also to density fluctuations (the cosmic web) and to the 
IGM temperature field (which affects the recombination 
rate and, hence, the He II fraction within each absorber). 

To account for density fluctuations in the He II forest, 
Fardal et al. (1998) considered variations in the counts 
of H I Lya forest lines together with the transformation 
-^Heii = ijNm- They found that the fractional standard 
deviation of the effective optical depth of the forest over 
a segment of length (in redshift space) Az is 3 
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where 7/ = 50?75o is the ratio of column densities for He II 
and H I and £hcII parameterizes the line width of the 
He II lines relative to H I: it is unity for turbulent broad- 
ening (or that provided by the Hubble flow) and one-half 
for thermal broadening. Theoretical models predict that 
the former is most relevant for the forest, but the mea- 
surements may be better fit by the latter (Fechner & 

3 Here we use the forest-dominated case of Fardal et al. (1998), 
rather than continuous IGM absorption. We also use their model 
A2 for the distribution of H I absorbers, which remains a reasonable 
fit to the low-column density absorbers relevant for the He II forest. 
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Reimers 2007b). Observations now show that (rj) w 50- 
80 (Heap et al. 2000; Shull et al. 2004; Zheng et al. 
2004b; Fechner et al. 2006; Fechncr & Reimers 2007a) 
with wide fluctuations between individual absorbing sys- 
tems (that themselves depend upon the ionizing back- 
ground and so are already incorporated in our model; 
Furlanetto 2009a). Fortunately, the extremely weak de- 
pendence on £hcII and 7750 implies that their uncertainties 
will not significantly impact our calculations. 

Temperature fluctuations are substantial at the end of 
He II reionization (Gleser et al. 2005; Furlanetto & Oh 
2008a; McQuinn ct al. 2009), and they will also affect 
the forest. However, we expect their effects to be signif- 
icantly smaller even than those due to density fluctua- 
tions. We have z H eii oc A 2 a(T)/r cx A 2 T-°- 7 /T, where 
A is the density in units of the mean and a(T) is the re- 
combination rate. The variance in T is significantly less 
than an order of magnitude (and hence, smaller than or 
comparable to that in A and T), even at the height of 
He II reionization, and is further damped by the smaller 
exponent. Thus that factor is much less important than 
the others, and we ignore it. 

Finally, our simple Monte Carlo model also ignores the 
additional fluctuations sourced by the intricacies of radia- 
tive transfer through the IGM (for example, shadowing 
and spectral filtering by dense regions), as well as fluctu- 
ations in the absorber population. For the latter, Zuo & 
Phinney (1993) showed that Poisson fluctuations in the 
absorber counts do not affect the correlation function of 
the ionizing background, so they should not affect our 
simpler statistics either. The details of radiative trans- 
fer are best addressed with more sophisticated numerical 
simulations (Maselli & Fcrrara 2005; Tittley & Meiksin 
2007; Paschos et al. 2007). We do not expect them to 
change our conclusions substantially, because we aver- 
age over large spatial regions and so are not subject to 
the large variations that radiative transfer can induce on 
small scales. 

2.4. The Attenuation Length 

With the quasar luminosity function now well- 
measured (especially at the bright end), the key input 
determining the overall amplitude of the fluctuations in 
our model is the attenuation length of ionizing photons. 
At z = 3, estimates in the literature range from <~ 30- 
60 Mpc at the ionization edge (Bolton et al. 2006; Furlan- 
etto & Oh 2008b; Faucher-Giguere et al. 2008), increas- 
ing by a factor of ~ 50% when averaged over the entire 
frequency spectrum. Most studies take Ro oc (1 + z)~ 3 
in the post-reionization limit, which we also assume as 
a fiducial model; this produces a reasonable match to 
the z < 2.7 He II Lya forest data (Dixon & Furlanetto 
2009). Of course, the overall level of fluctuations in- 
creases as the attenuation length decreases (Zuo 1992a; 
Fardal & Shull 1993; Meiksin & White 2003), although 
the effect is not dramatic because of the wide range of 
intrinsic quasar luminosities (Meiksin 2007; Furlanetto 
2009a) . We will consider two fiducial attenuation lengths 
below, although we will not account for the much larger 
attenuation lengths of high-energy photons. These will 
smooth out the background and so only strengthen our 
argument that the "null hypothesis" of a smoothly evolv- 
ing attenuation length is inconsistent with the observed 
fluctuation amplitude at z ~ 3. 



3. A TOY MODEL FOR REIONIZATION 

So far, we have assumed that He II is highly-ionized 
everywhere, so that the only opacity source for high- 
energy photons is the (uniform) attenuation represented 
by i?o- During reionization, however, large swathes of 
the IGM contain primarily He II, absorbing nearly all 
photons above the ionization edge. Such regions will, 
of course, appear as complete absorption troughs in the 
spectrum, naturally increasing the expected fluctuation 
amplitude. 

Unfortunately, a detailed model of the morphology of 
He II and He III regions requires a full numerical simu- 
lation of the reionization process, as in McQuinn et al. 
(2009). Rather than attempt such a calculation, we will 
instead use a simple toy model to estimate the impact of 
reionization on the He II Lya forest variations; we will, 
therefore, forego precise quantitative constraints during 
this era and, instead, use this model to demonstrate that 
reionization can plausibly explain many aspects of the 
observations. 

Specifically, we first choose the fraction Xdark of pixels 
that lie outside of regions illuminated by quasars. This 
quantity is nearly, but not quite, the He II fraction XHeii: 
^dark > ^Heii because some regions that were ionized 
by earlier generations of quasars will no longer have a 
source illuminating them. Without such a nearby source, 
recombinations are sufficiently rapid to quickly render 
t cS > 1 (Furlanetto et al. 2008; Furlanetto 2009a). 

We then assume that these opaque regions are precisely 
those with the smallest T in the post-reionization limit; 
in other words, they are the regions least illuminated 
by existing quasars. We identify the threshold Tdark for 
which the probability that a random pixel has T < Tdark 
as rrdark- We then set T = for all pixels below this 
threshold. We emphasize that this is an overly simplified 
prescription, because it does not account for the pos- 
sibility that a He II region may be relatively close to 
an existing quasar. It is accurate in the limit in which 
quasars are extremely long-lived. 

Despite this rather extreme simplification, the model 
does have the nice property that brightly-illuminated re- 
gions near existing quasars are preserved (unlike if, for 
example, we laid down He II regions randomly), which is 
qualitatively consistent with the evolving distribution of 
r throughout reionization (Furlanetto 2009a). Moreover, 
given that previously ionized regions will quickly become 
opaque due to recombinations (Furlanetto et al. 2008), 
this simple model does provide the correct qualitative 
distribution. It also does not require us to prescribe the 
spatial sizes of ionized regions. 

Figure 1 shows some example lines of sight generated 
with this toy model if rrdark = 0-2- The three original 
(post-reionization) lines of sight are shown in the bot- 
tom panel, to make clear the effects of our toy model. 
One interesting aspect of the prescription is immediately 
obvious: the dark regions span a wide range of scales, 
from <~ 10 Mpc to several hundred Mpc (from r ~ 250- 
450 Mpc along the dashed curve, for example). Thus, at 
least in this toy model, both the early and late stages 
of reionization will enhance the observed fluctuations in 
r e ff, even when we smooth the data on very large scales. 

4. EFFECTIVE OPTICAL DEPTH DISTRIBUTIONS 
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^e,f/< T e,,> T .,f/< T e,f> 

Fig. 2. — Left: Probability distributions of r e jf averaged over 100 comoving Mpc segments in our simulations, scaled to an imposed mean 
value. The three panels take different points during and after He II reionization; the different curves normalize (r e f{) differently, as labeled 
in the bottom panel. All are measured from our Monte Carlo model at z = 3, assuming Rq = 30 Mpc. Right: Same, but showing the 
cumulative probability distributions. 



The left panels of Figure 2 show the r e s probability dis- 
tributions for several simple cases. Each takes segments 
of 100 Mpc (comoving), z = 3, Rq = 30 Mpc, and uses 
our full Monte Carlo sample (with 3x1 Gpc 3 volumes, 
each with 100 lines of sight 1 Gpc long). In each panel, 
we allow three different normalizations for (r e ff). The 
right panel shows the cumulative distribution functions. 
We ignore density fluctuations in these models. 

The post-reionization distributions (in the bottom left 
panel) contain substantial, but not overwhelming, varia- 
tions, with a full width at half maximum comparable to 
the mean value. In agreement with Fardal et al. (1998), 
we find that the scaled distributions are relatively insen- 
sitive to (T er r); our variance is also comparable to theirs, 
though our different choice for i?o has some effect. A 
larger mean opacity introduces more data points with 
moderately large T e g, because it is easier for a long IGM 
segment to lie entirely at low T. However, the cumu- 
lative distributions make clear that none of these cases 
have regions with r e s/ (T e s) > 2.5 at a very high level of 
confidence (i.e., none appear in our 3,000 samples). 

Figure 3 quantifies the post-reionization variations in 
a slightly simpler way. Here, we show the medians of 
the distributions, as well as the limits within which 68% 
and 95% of the simulated segments lie. First, note that 
the median is generally a good estimate of the average, 
although the skewness at larger (r e g) does introduce a 
small bias. The distribution remains relatively invari- 
ant toward small optical depths throughout the range in 
r ff. Physically, strong transmission occurs near bright 
quasars, where the strong peaks in T are independent of 
the average transmission. On the other hand, the high- 
r e ff tail does broaden substantially at larger mean optical 
depths. 

Figure 3 also shows the impact of Rq on the distribu- 
tions; the thick and thin curves assume Rq = 30 and 50 
Mpc, respectively. Overall, a larger attenuation length 



means more sources contribute to T at each point in the 
IGM, so the stochastic fluctuations are smaller - and, 
indeed, the distribution of t c s narrows somewhat as Rq 
increases. However, the effect is not symmetric: it is 
quite modest at r e g < (r e s), because such segments are 
typically dominated by one or more luminous quasars, 
in whose immediate environs attenuation is unimportant. 
On the other hand, the variations with Rq are fairly large 
at T c ff > (r e ff). They are especially significant at the far 
end of this tail (the upper dotted curves), because it be- 
comes very difficult to find long regions with extremely 
small values of T as the volume sampled by each IGM 
point increases. 

The other panels in Figure 2 show the distributions 
from our toy model of He III reionization (for Xdark = 0.2 
and 0.5 for the middle and top panels, respectively; in the 
latter case, we show only (r c ff) = 2 and 3, since a lower 
mean optical depth so early in reionization seems implau- 
sible). Even in these cases, the bulk of the distribution 
does not change dramatically from the post-reionization 
limit. However, a long tail toward high opacity clearly 
does develop; the cumulative distribution shows that 
~ 5% and 15% of the 100 Mpc skewers are completely 
dark for £dark =0.2 and 0.5, respectively. 4 

Note that the observed fraction of dark segments is 
much smaller than Xdark itself - because even a single 
bright quasar near the 100 Mpc skewer suffices to pro- 
duce measurable transmission, if only over a short region 
- but, even over long segments, a non-negligible fraction 
of completely dark segments remains at the end of reion- 
ization (see Fig. 1). Thus, at least in our simple models, 
reionization clearly can increase the measured variance, 
primarily by generating a tail toward extremely strong 
absorption, even after averaging the segments over large 

4 Actually, these regions will still have a weak ionizing back- 
ground from X-ray photons, but this radiation field is not strong 
enough to render them transparent (Furlanetto 2009a). 
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Fig. 3. — Scaled distributions of measured r e ff values as a func- 
tion of the true average opacity in the post-reionization limit. The 
distributions average over 100 Mpc segments at z = 3 and assume 
Ro = 30 and 50 Mpc (thick and thin curves, respectively). The 
solid curves show the medians of the distribution; the dashed and 
dotted curves show the limits within which 68% and 95% of the 
simulated measurements lie. 



Fig. 4. — Compiled He II Lya forest effective optical depth 
measurements from the literature (see Dixon &: Furlanetto 2009 
for references) . The two curves are predictions of the optical depth 
assuming a smoothly evolving mean free path and the Hopkins 
et al. (2007) quasar luminosity function, normalized to the observed 
(r e ff) at z = 2.45. The solid and dotted curves assume Ro = 30 or 
60 Mpc at 2 = 3, respectively. 



path lengths. 

In addition, when (r c g) is relatively small, a more 
pronounced asymmetry does develop in the distribution, 
with the median falling somewhat farther below the av- 
erage value (see, e.g., the solid curve in the middle panel 
on the left) . These cases have an exceptionally high con- 
trast between the dark regions and those around lumi- 
nous quasars (which provide almost all of the transmis- 
sion). 

So far, we have ignored fluctuations from the density 
structure (manifested through the line structure of the 
Lya forest); as shown by Fardal et al. (1998), they are 
typically ~ 0.1 (r c ff) for segments of length 100 Mpc (or 
Az w 0.1; see eq. 7). This is indeed a few times smaller 
than those from the radiation field, unless Rq is many 
times larger than our fiducial values. 

5. COMPARISON TO DATA 

Although (to date) only five lines of sight have detailed 
He II Lya forest spectra, the data still contain substan- 
tial fluctuations to which we can compare our model. 
Figure 4 shows all of the (averaged) r c g measurements 
at z < 3.2 from the literature, as compiled by Dixon & 
Furlanetto (2009). The data come from HE 2347-4342 
(Kriss et al. 2001; Shull et al. 2004; Zheng et al. 2004b), 
HS 1700+64 (Fechner et al. 2006), Q0302-003 (Hogan 
et al. 1997; Heap et al. 2000; Worseck & Wisotzki 2006), 
HS 1157+314 (Reimers et al. 2005), and PKS 1935-692 
(Anderson et al. 1999). 5 We have binned the measure- 
ments over Az = 0.1; within each bin, we show one data 
point for each line of sight, as well as the average over 
all the lines of sight (solid squares). 

The figure also shows two predictions for the evolu- 
tion of T e ff assuming that helium reionization occurred 

5 We do not include the quasar discovered by Zheng et al. 
(2004a), because it provides only a lower limit to r c ff at z ~ 3.5. 



at z > 3.2 (see Dixon & Furlanetto 2009 for details). 
Both curves take as inputs the quasar luminosity func- 
tion from Hopkins et al. (2007) and a prescription for 
the attenuation length; we assume Rq oc (1 + z)~ 3 with 
i?o = 30 and 60 Mpc at z = 3 for the solid and dotted 
curves, respectively. These two values span the range of 
estimates in the literature for photons at the ionization 
edge (e.g., Bolton et al. 2006; Furlanetto & Oh 2008b; 
Faucher-Giguere et al. 2008). 

With these input parameters, Dixon & Furlanetto 
(2009) combined the fluctuating radiation field models of 
Furlanetto (2009a) and the density field fit from Miralda- 
Escude et al. (2000) to generate (r e s) via the fluctuating 
Gunn-Peterson approximation. Because this approxima- 
tion does not fully account for the line structure of the 
Lya forest, we have normalized both curves to the data 
point at z = 2.45, which we regard as relatively robust. 
This method allows us to compare the predicted opti- 
cal depth evolution to the observations, irrespective of 
any concerns about the normalization of the emissivity, 
the recombination rate, or the mean free path. Because 
of this normalization, the two models yield very similar 
overall evolution for (r ff), although of course the model 
with the larger attenuation length implies smaller over- 
all fluctuations in the radiation field (and hence smaller 
variations in T e s). 

There are several regimes that contain interesting dis- 
crepancies with our smooth model, as described in Dixon 
& Furlanetto (2009), and we will examine each in detail 
below. First, at z < 2.25, the data from HE 2347-4342 
imply that r is ~ 2 times smaller than our simple ex- 
trapolation from higher redshifts would imply. Second, 
at z ~ 2.7-2.9, the line of sight to HE 2347-4342 shows 
substantial fluctuations between large and small trans- 
missions on > 10 Mpc scales. Finally, at z > 2.9, the line 
of sight to Q0302-003 shows a substantial increase in r e g 
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Table 1. Probability Estimates for Ro(z = 3) = 30 Mpc 



<*> 


Ar (Mpc) 


fcff (obs) 


(r cff ) (model) 


P(lEdark = 0) 


P(z d ark = 0.2) 


P(Zdark = 0.5) 


2.05 
2.15 
2.25 


146 
140 
134 


0.4991°;°^ 
0.589±°°il a 
0.837t°;° 3 34 a 


0.455 
0.520 
0.605 


30% 
26% 
3% 


_b 
_b 
_b 


_b 
_b 
_b 


2.76 
2.83 


35 
20 


o 7 -+U.V«c,d 
°- ,o -0.43 
q 1 o0.58 c 
• 3 - 1 °-0.37 


1.28 
1.50 


< 0.01% 
1% 


10% 

16% 


35% 
33% 


2.93 
3.07 
3.15 


61 
58 
94 


4 - J '-0.51 

z - oo -0.09 
4 74 +0.46c 
4 - ,4 -0.31 


1.72 
1.90 
2.23 


< 0.03% 

34% 
2% 


7% 
33% 
6% 


25% 
45% 
17% 



a From HE 2347-4342 (Zheng ct al. 2004b). 

b We assume that helium reionization is complete by z = 2.3 so do not report these values. 
c From HE 2347-4342 (Shull ct al. 2004). 

d Note that Kriss et al. (2001) report a lower limit in this region of T e g > 2.75 with the same data. 

°From Q0302-003. Measurements are from Worseck & Wisotzki (2006), using data from Heap et al. (2000). 



over the smooth models. Thus we will focus on these two 
lines of sight (amongst the best observed) as particular 
instances of interesting fluctuations. Our primary goal 
is to determine whether the data can be explained by 
including cosmic variance around the smoothly evolving 
model shown here. If not, they imply substantial changes 
in the input parameters - most likely He II reionization. 

Throughout the following, we will ignore variations in 
the underlying density field; equation (7) shows that such 
deviations are modest and do not substantially affect our 
conclusions. 

5.1. Fluctuations at z < 2.3 

We first wish to assess the significance of deviations 
from the "smooth" T evolution models at low redshift. 
Near z ~ 2, Dixon & Furlanetto (2009) found that 
i?o oc (l + z)~ 3 overpredicted the apparent ionizing back- 
ground - although only one line of sight has provided 
measurements in this regime (HE 2347-4342). The de- 
viations are much larger than the statistical errors from 
the observations; although in absolute terms, they only 
require a factor of two decrease in T. 

We list the measured optical depths from HE 2347- 
4342 (Zheng ct al. 2004b) and the predictions from our 
two models at z < 2.3 in Tables 1 and 2. The first four 
columns in each table give the central redshift of the 
skewers, the comoving width of the skewers, the observed 
mean transmission (and la statistical errors), and the 
predicted effective optical depth, respectively. 

Before evaluating whether these points are truly dis- 
crepant with the models, we do note two potential sys- 
tematic problems with the data. First, this low-z regime 
corresponds to small wavelengths, so the continuum cor- 
rection is most uncertain here. Second, this regime has 
been corrected for Ly/3 and higher-order line contamina- 
tion using the corresponding Lya forest, but that cor- 
rection may not be perfect given the noisy He II data. 
Neither of these factors are included in the error esti- 
mates, so they already decrease the significance of the 
apparent discrepancy. 

To estimate the significance of these deviations, we use 
our Monte Carlo model to generate 100 lines of sight in 
each of three different realizations of the quasar distri- 
bution. Each box is ten times larger than the Az = 0.1 
segments used in the data binning. We thus have 3000 



total segments for each comparison. We report the frac- 
tion of these model segments that have r G ff greater than 
or equal to the observed values in the fifth columns of 
Tables 1 and 2. 

Clearly, the relatively small deviations at z < 2.2 are 
not at all surprising in either model. However, the de- 
viation at z — 2.25 is significant, occurring in at most 
a few percent of the lines of sight. Such a discrepancy 
could more easily be accommodated if we allow fully dark 
regions (i.e., demand that He II reionization is still ongo- 
ing), but we do not regard that as anywhere near likely. 

One obvious solution is that the model (r e ff) (z) falls 
off somewhat too steeply with redshift compared to the 
data. 6 Only a modest change would be necessary; in- 
creasing the model prediction by A (r e ff) = 0.05 increases 
the expected probability to ~ 15% in the Rq = 30 Mpc 
model. Indeed, Fechner et al. (2006) find that the HS 
1700+64 data are consistent with a smoothly evolving 
power law somewhat flatter than our curves, which would 
be completely consistent with the lower-redshift points 
(but make the disagreement at higher redshifts even more 
apparent). 

Alternatively, density fluctuations may help reconcile 
the data. Equation (7) suggests that this mechanism 
provides At c s ~ 0.05. This kind of uncertainty can ul- 
timately be calibrated by comparing to the H I Lya for- 
est. However, in the case of the z ~ 2.25 region toward 
HE2347-4342, a visual inspection shows if anything more 
H I transmission than average (Zheng et al. 2004b) . 

One final worry is that all three data points at z < 2.3 
show a systematic increase in r e s relative to the simple 
model, whereas one might expect fluctuations on both 
sides of the mean. However, Figure 1 shows that it is 
not particularly uncommon to find regions several hun- 
dred Mpc long over which T remains small. To examine 
the likelihood of this, we consider estimates of the trans- 
mission over segments of length 420 Mpc at z = 2.15. We 
find that 7% (2%) of the model regions have r c g above 
the observed value if Rq = 30 (60) Mpc. These arc some- 
what more likely than the z = 2.25 segment on its own, 
so such large scale deviations are not particularly sur- 

6 As described in detail by Dixon & Furlanetto (2009), another 
possibility is an evolving normalization factor in the fluctuating 
Gunn-Peterson approximation used there. 
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Table 2. Probability Estimates for Ro(z = 3) = 60 Mpc 



<*> 


Ar (Mpc) 


fcff (obs) 


(r cff ) (model) 


P(lEdark = 0) 


P(z d ark = 0.2) 


P(Zdark = 0.5) 


2.05 
2.15 
2.25 


146 
140 
134 


0.499t°;°^ 
0.589±°°il a 
0.837t°;° 3 34 a 


0.473 
0.538 
0.617 


30% 
26% 
0.03% 


_b 
_b 
_b 


_b 
_b 
_b 


2.76 
2.83 


35 
20 


o 7 -+U.V«c,d 
°- ,o -0.43 
q 1 o0.58 c 
• 3 - 1 °-0.37 


1.25 
1.45 


< 0.01% 

< 0.01% 


10% 

20% 


36% 
35% 


2.93 
3.07 
3.15 


61 
58 
94 


4 - J '-0.51 

z - oo -0.09 
4 74 +0.46 e 
4 - ,4 -0.31 


1.66 
1.96 
2.17 


< 0.03% 
18% 

< 0.03% 


9% 
22% 
6% 


28% 
45% 
20% 



a From HE 2347-4342. 

b We assume that helium reionization is complete by z = 2.3 so do not report these values. 
c From HE 2347-4342 (Shull ct al. 2004). 

d Note that Kriss et al. (2001) report a lower limit in this region of T e g > 2.75 with the same data. 

°From Q0302-003. Measurements are from Worseck & Wisotzki (2006), using data from Heap et al. (2000). 



prising. 

5.2. Fluctuations at z ~ 2.8 

The second interesting region is at z <~ 2.8, where the 
HE 2347-4342 line of sight shows substantial fluctuations 
(both on its own and relative to the other lines of sight). 
We first focus on the segment at z = 2.76, which spans 
~ 35 Mpc and is nearly completely opaque (as measured 
most recently by Shull et al. 2004). Here, we again use 
three independent boxes, 100 lines of sight in each box, 
and 30 short segments along each line of sight (for a total 
of 9,000 samples). Tables 1 and 2 show that neither of our 
post-reionization models can accommodate this point at 
a confidence level of < 0.01%. Even if wc use the reported 
3cr lower limit r e s > 2.78, the probability is still less than 
0.01% in cither model. 

On the other hand, if we allow for dark regions from in- 
complete helium reionization, the model probability be- 
comes reasonably large: > 10% (35%) for rzidark = 0.2 
(0.5). Thus, discounting any possible systematic errors 
in the data, this single point is extremely suggestive that 
helium reionization has not yet completed. 

Physically, this conclusion is so powerful, because the 
post-reionization models have a sharply defined mini- 
mum in T throughout the Universe, generated by the 
accumulated background of quasars at large impact pa- 
rameters (within two or three attenuation lengths) . This 
level is > 0.31", so it is very difficult to accommodate 
large optical depths, even over relatively modest path 
lengths. (On the other hand, excursions to small optical 
depths are easy to accommodate by passing near a bright 
quasar.) 

One possible explanation is a systematic error in the 
measurement. The relevant portion of the FUSE spec- 
trum is relatively close to an airglow emission line at 1136 
A (corresponding to z = 2.73), but overall this region is 
no more or less difficult to analyze than others in the 
spectrum (J. Tumlinson, private communication). Kriss 
et al. (2001) examined the same data set and placed a 
la lower limit of r e ff > 2. 6, 7 perhaps indicating the level 
of systematic uncertainties on the background subtrac- 
tion, etc. Nevertheless, the probability of finding such a 

7 Note that this is averaged over a slightly different wavelength 
interval. 



segment in our post-reionization models is still < 0.1%. 

Another possibility is that (r e ff) evolves more steeply 
at high redshifts, so that the mean value is actually closer 
to the measured value. Figure 4 shows that the measured 
average T e g is also somewhat higher than our model. 
However, if we assume that (t c s) = 1-75 at z — 2.76, 
somewhat above the measured average, we find that the 
probabilities for r e s > 3.75 are < 0.01%, although the 
probabilities for r e s > 2.75 are more reasonable (<~ 10% 
and 1% for the two mean free paths). Thus a combi- 
nation of the Kriss et al. (2001) analysis and an under- 
estimate in the mean from our model would provide a 
plausible scenario. 

However, as argued in Dixon & Furlanetto (2009), it 
would be difficult to explain such a sudden upturn with- 
out appealing to He II reionization. The only possible 
causes are (1) rapid evolution in the ionizing emissiv- 
ity or (2) rapid evolution in the attenuation length. The 
first is ruled out by observations of the quasar luminosity 
function (e.g., Hopkins et al. 2007). The second would 
probably be associated with the continued clearing of 
IGM He II during reionization, unless the modest in- 
crease in r associated with reionization filters through 
the moderately dense He II absorbers over a long time 
scale. 

Of course, effects we have ignored in our model can 
also moderate the discrepancy somewhat. For example, 
including density fluctuations (which, in our model, man- 
ifest themselves through the number of Lya forest ab- 
sorbers in a segment) increases the variations slightly, 
though not sufficiently to reconcile this measurement 
with the model. More promising is radiative transfer: 
one can imagine that shadowing, where a dense system 
blocks the light from a nearby source, could produce a 
regions with small T and hence little transmission. Still, 
because the lower limit on T is provided by the collective 
action of many distant quasars and because the segment 
is relatively long (~ 35 Mpc) , it seems difficult to arrange 
strong shadowing over the entire interval. 

The same line of sight also contains strong fluctuations 
at z = 2.8-2.88 toward both large and small transmission 
(Kriss ct al. 2001; Shull et al. 2004; Zheng et al. 2004b), 
with typical scales > 10 Mpc. 8 These features suggest 

8 These features are not shown in Figure 4 because they are 
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similar conclusions, albeit at somewhat lower confidence. 
In the table, we focus on a single measurement from Shull 
et al. (2004) at z « 2.83 that has a relatively large path 
length; the other strong fluctuations (toward both large 
and small transmission) are more easily accommodated 
because they have shorter path lengths. 

5.3. Fluctuations at z > 2.9 

We next consider the z > 2.9 range probed by Q0302- 
003. In Dixon & Furlanetto (2009), we found that these 
data imply a sharp break in (r e fj) at z < 2.9, with 
substantially larger values at higher redshift. For ref- 
erence, we provide the binned t c s measurements from 
Worseck & Wisotzki (2006), who re-analyzed the Heap 
et al. (2000) observations, in Tables 1 and 2. At z = 3 
(corresponding to 1216 A), geocoronal H I Lya contami- 
nates the spectrum; we have excised the range z = 2.96- 
3.04 to conservatively account for this. At z — 2.8-2.9, 
the opacity along this line of sight is consistent with 
other measurements. Above this range, the absorption 
increases sharply, then decreases (thanks largely to a 
nearby quasar discovered after the spectrum was taken; 
Jakobsen et al. 2003), and resumes increasing. 

The last sections of Tables 1 and 2 show the likeli- 
hood of finding the measured values in the context of 
our smoothly evolving models. 9 Clearly, the point at 
z = 2.93 is strongly discrepant with the post-reionization 
models, but it can be relatively easily accommodated if 
a^dark ^ 0.2. The same is true for the z = 3.15 point. On 
the other hand, at z = 3.07, the substantial transmission 
spike from the nearby quasar makes the data entirely 
consistent with the smoothly evolving models. 

Of course, one can easily accommodate the measure- 
ments if our model underestimates (r e ff) at z > 2.9 by 
a significant factor, but as we have argued above, this 
itself would seem to require ongoing He II rcionization. 
In the end, reconciling the high- and low-z behavior into 
smooth evolution appears difficult, requiring at least one 
model parameter to vary significantly more rapidly than 
typically assumed. 

These conclusions are perhaps not surprising, given the 
other evidence indicating that He II reionization occurs 
at z > 3. However, the fact that two high-r c ff regions 
appear (out of four possible), when the probability of 
each is < 10% at a^ark ~ 0.8, perhaps also indicates 
that even z — 3 is deeper in the reionization era than 
previously supposed. 

6. DISCUSSION 

We have computed the fluctuations expected in mea- 
surements of the average transmission of the He II Lya 
forest due to large-scale variations in the radiation field. 
We used a simple Monte Carlo model of the quasar dis- 
tribution to simulate the radiation field. We then mea- 
sured r e ff averaged over long (~ 30-100 Mpc) segments 
along random lines of sight and computed its probability 

small-scale, and the average opacity over the entire z = 2.8-2.88 
segment toward HE2347-4342 is consistent with just a modest up- 
turn from lower redshifts and is also consistent with other lines of 
sight at the same redshifts. Nevertheless, the fluctuations are quite 
dramatic; see e.g. Figure 3 of Dixon & Furlanetto (2009). 

9 We do not list the z = 2.8—2.9 region here because it is com- 
pletely consistent with the smooth models and with the other lines 
of sight. 



distribution. For a fully-ionized Universe, the radiation 
field typically produces a full-width at half-maximum for 
the r e g distribution of order its mean value (Fardal et al. 
1998). 

We also constructed a toy model of reionization, where 
we demanded that a fraction Xdark of pixels have zero 
transmission. We chose these dark pixels by identifying 
the regions that would have the smallest transmission in 
a post-reionization Universe. Such a model is reason- 
able, because only active quasars can illuminate a region 
strongly enough to render it transparent, even if fossil, 
mostly-ionized regions fill much of the IGM (Furlanetto 
et al. 2008). This method produces long dark stretches 
(with sizes > 100 Mpc) in the empty voids between 
quasars, even late in reionization. This procedure broad- 
ens the r e ff distribution by producing a large tail toward 
very high optical depths (and, to compensate, shifting 
the median transmission to smaller values) . 

We then compared our model to the observed fluctu- 
ations in the existing sample of He II Lya forest lines 
of sight, as compiled by Dixon & Furlanetto (2009). In 
particular, we asked whether a post-reionization model 
with a smoothly evolving attenuation length (and hence 
no change in the ionization state of IGM helium) could 
account for the observed data. 

At lower redshifts (z <~ 2.25) such a model is rea- 
sonable, even though the data show somewhat stronger 
absorption than expected. The observed transmission 
is unlikely according to our fiducial model, but modest 
changes to the normalization or slope of the underlying 
(r) evolution can easily accommodate it. Alternatively, 
the underlying IGM density fluctuations may also ac- 
count for the discrepancy (Fardal et al. 1998). 

This points to one important (though unsurprising) 
lesson from our analysis: even with the modest signal- 
to-noise already available in He II Lya forest spectra, the 
statistical errors on t c s arc significantly smaller than the 
natural level of cosmic variance. Thus, it is crucial to 
consider the latter when interpreting the existing data. 

At higher redshifts, the data become much more diffi- 
cult to explain. Even though there are only five lines of 
sight with existing measurements, we found strong evi- 
dence that a smoothly-evolving post-reionization model 
is inconsistent with the observations. At z ~ 2.7-2.9, 
the measurements fluctuate strongly toward high opti- 
cal depths on ~ 20-35 Mpc scales. Such regions of deep 
absorption cannot fit into a post-reionization model (at 
> 99% confidence) , because the large number of quasars 
far from the line of sight already provide enough radia- 
tion to keep the ionized fraction reasonably large. Thus, 
the deviation cannot simply be explained as cosmic vari- 
ance from a single line of sight, unless it is truly patho- 
logical. 

To reconcile the observations and model, we need to 
increase the predicted (r e ff) significantly, use relatively 
small mean free paths, and take a conservative view of 
the data (i.e., accepting the more conservative Kriss et al. 
2001 lower limits). However, in this case there still must 
be a break in (r e ff) (z) at z ~ 2.7, which we have previ- 
ously argued would likely itself signal He II reionization 
(Dixon & Furlanetto 2009). 

At z > 2.9, the only well-observed line of sight (Q0302- 
003) has r e ff > 4 in two separate ~ 60-100 Mpc intervals. 
This opacity is much larger than our smoothly-evolving 
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model predicts, and even when post-reionization radia- 
tion fluctuations and density fluctuations are included, 
it remains extremely improbable. Thus, this line of sight 
on its own also argues strongly for either ongoing reion- 
ization or a rapid increase in (r e ff) at z <~ 2.9, likely also 
signaling the end of reionization. 

Our model, therefore, supports the conclusion of Dixon 
& Furlanetto (2009), who examined the evolving average 
effective optical depth, that the data are best explained 
by models with reionization completing at z < 2.9. This 
is somewhat later than expected from measurements of 
the IGM temperature (Schaye et al. 2000; Ricotti et al. 
2000; Lidz et al. 2009). It is more consistent with 
hints that the shape of the ionizing background spec- 
trum changes around that time (Songaila 1998, 2005), 
although that argument has remained controversial (Kim 
et al. 2002; Aguirre et al. 2004). Late reionization is rea- 
sonably consistent with estimates of the total ionizing 
emissivity of high- z quasars, which produce a couple of 
ionizing photons per helium atom by this point (e.g., 
Furlanetto & Oh 2008b). 

We have focused on variations in the radiation back- 
ground. Fardal et al. (1998) showed that fluctuations 
induced by large-scale IGM density variations are mod- 
est. While they will slightly modify the values in Tables 1 
and 2, they will not qualitatively modify our conclusions 
about the high-r c ff measurements. They can, however, 
help to reduce the tension at z < 2.3, where the absolute 
deviations are relatively small. 

Our toy model of reionization can obviously be much 
improved with realistic simulations of reionization, such 
as those performed by McQuinn et al. (2009) - although 
the large volume required to accumulate enough sam- 
ples presents a challenge. The probabilities generated 
from our simple model only provide a plausibility argu- 
ment that this event, where large coherent regions with 
weak radiation fields should be relatively common, can 
more easily accommodate the observed T c g distributions. 
Qualitatively, we expect similar behavior in more sophis- 
ticated models (see, for example, the visualizations in 
Figs. 4-7 of McQuinn et al. 2009). We therefore have 
not tried to estimate the evolution of the He III fraction 
with redshift but instead to argue that the data suggest 
some sort of change in the ionized fraction at z < 2.9. 

Our post-reionization model does require some caveats. 
First, it does not include the details of radiative trans- 
fer, such as shadowing and spectral filtering. These ef- 
fects can clearly cause substantial fluctuations across the 
absorber population (Maselli & Ferrara 2005; Tittley & 
Meiksin 2007). However, by averaging over large seg- 
ments of each spectrum (> 30 Mpc across), we should 
also average over these details. For example, it is diffi- 
cult to imagine an absorber large enough to shadow an 
entire 30 Mpc region from the many sources that illumi- 
nate it within distances of a few attenuation lengths. 

Second, we have ignored the deterministic clustering 
of quasars (Shen et al. 2007). Furlanetto (2009a) showed 
that this will be relatively unimportant to the distri- 
bution of r, because the fractional variance in quasar 
emissivity within one attenuation volume is much larger 
than that produced by deterministic clustering (thanks 
to the large Poisson fluctuations and the wide luminos- 
ity range of quasars). McQuinn et al. (2009) also showed 
that the distribution of He II during reionization is es- 



sentially driven by the Poisson fluctuations of the quasar 
population. However, quasar clustering may be some- 
what more important for our application, because clus- 
tering within the attenuation volumes cannot really be 
neglected (Mesinger & Furlanetto 2009) and because the 
spatial coherence is important for averaging over long 
skewers. However, it should still not be a dominant ef- 
fect, and we expect that our qualitative conclusions will 
hold. 

We have also simplified the problem by taking a sin- 
gle attenuation length for all photons, whereas in reality 
it increases with frequency. However, the pair of val- 
ues we have used spans a factor of two in length, both 
yielding similar conclusions, so we do not expect this to 
make a large difference. If anything, including higher- 
energy photons will strengthen evidence for late reion- 
ization, because larger attenuation lengths decrease the 
fluctuations. 

It may seem surprising that we have drawn such strong 
conclusions, given the very limited number of lines of 
sight so far available. The most conservative conclusion 
one could draw is that the data are inconsistent with a 
smoothly evolving ionizing background over the entire 
range z ~ 2-3.2. Either (1) there must be a sudden in- 
crease in the slope of (r e ff) (z) or (2) extra fluctuations 
from reionization are required. These conclusions are 
robust because they rely on the difficulty of producing 
a region of any substantial length with r <C T in the 
post-reionization Universe. Thus, so long as an interval 
is large enough to include a substantial number of ab- 
sorbers (and so average over the density fluctuations), 
the detection of even a single region with r e ff 3> (T e ff) 
provides strong evidence for reionization. 

In any case, the next several years should see a sub- 
stantial increase in the number of lines of sight, with 
dozens of promising targets recently detected (Zheng 
et al. 2004a, 2008; Syphers et al. 2009a,b). The instal- 
lation of the Cosmic Origins Spectrograph and repair of 
the Space Telescope Imaging Spectrograph on the Hub- 
ble Space Telescope (HST) have provided two powerful 
instruments for such exploration; intriguingly, these in- 
struments are best suited for studying the IGM at z > 3, 
precisely where the fluctuations should be most inter- 
esting. Syphers et al. (2009a,b) have detected 22 prime 
candidates for He II Lya forest spectra. Nearly all lie 
in the extremely interesting redshift range z ~ 3-4, in 
the heart of the reionization era. A visual inspection of 
the modest-resolution spectra available from their sur- 
veys (taken with the ACS grism on HST) shows many 
interesting features both within individual spectra and 
across the entire population. 

Finally, our model shows that even modest-resolution 
and relatively low signal-to-noise investigations (such as 
the ACS grism sample) can make substantial progress in 
understanding He II reionization. We have drawn strong 
conclusions from the absorption averaged over fairly large 
intervals (Az ~ 0.03-0.1), which is (at least in principle) 
a resolution- independent measurement. Thus, this test 
can be useful even for quasars that are too faint for the 
detailed He II Lya forest measurements that have so far 
held the focus of the community. 
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